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INTERACTION OF STRUCTURE AND LIQUID 
IN THE SOUND SUPPRESSOR SYSTEM 

SUMMARY 


An analysis of the interaction of structure and liquid of the Saturn V 
sound suppressor system has been performed. For this reason the liquid re- 
sponse in a rectangular container to translatory excitation has been determined 
and was described by a mechanical analogon. It was found that more than eighty 
percent of the liquid mass in the water reservoir sloshes for the fundamental 
sloshing mode. The natural frequencies of the liquid are extremely low; the 
fundamental frequency is in the vicinity of 0. 2-0. 3 or 0. 4-0. 5 rad/sec depend- 
ing on the orientation of the excitation. An analysis of the interaction of the 
liquid and the sound suppressor structure revealed rather low coupled frequen- 
cies of 0. 22 and 0. 74 rad/sec in longitudinal direction (y-direction) and about 
0. 5 and 6. 7 rad/sec in the cross direction ( z- direction ) . It was found that an 
increase of stiffness does not pay off. Because the cross direction offers the 
largest structural area, the response of the system to some winds has been 
treated. Liquid and structural amplitudes are presented for a sinusoidal gust 
wind of 13 m/sec amplitude, for rectangular pulses of various durations as well 
as for exponentially decaying pulses of different decay magnitude. 


SECTION I. INTRODUCTION 


Through model tests employing live rocket engines it has been determined 
that the most reliable and effective means to suppress the high-intensity low- 
frequency noise generated by large space vehicle booster engines during test 
firings is the injection of large quantities of water near the engine nozzle exit. 

This reduces the kinetic energy of the high velocity gases. A reservoir of about 
five million gallon capacity is necessary to supply the water for the suppression 
of sound of the Saturn V rocket engines. The water must be supplied just below 
the elevation of the engine exit plane which is almost 100 feet above grade (Fig. i). 

The operation of the sound suppressor or some wind forces on the total 
system may cause destructive forces or uncontrollable conditions which could 
lead to the malfunction of the sound suppressor or even to the destruction of the 



total system. It is, therefore, necessary to investigate the motion of the liquid 
in the container and its coupling with the elastic support structure to reveal the 
integrity of the intended design and make recommendations for a final enhanced 
system. 

In the following, the response of the liquid in the container due to trans- 
latory excitation is determined and is described as a mechanical analogon. From 
this investigation it will be learned how many of the liquid vibration modes will 
have to be included in the overall structural analysis. After the equations of 
motion of the system have been derived, the coupled frequencies will be deter- 
mined for various stiffnesses of the support structure. Forced undamped and 
damped harmonic response of the system is studied and finally the response of 
the system due to various wind inputs is derived. 


SECTION II. RESPONSE OF LIQUID TO TRANSLATORY EXCITATION 


For the structural analysis of the sound suppressor system the response 
of the liquid in the reservoir due to some outside excitation has to be known. It 
is also important to know its natural frequencies. Because of the system design 
only translatory excitation of the liquid container is probable. For this reason, 
investigation is restricted to this type of excitation for a rectangular container 
of infinite width (Fig. 2). If the liquid is considered incompressible, irrota- 
tional and nonviscous, the velocity of the liquid can be described by the gradient 
of a velocity potential, which because of the continuity equation must be a solution 
of the Laplace equation 


V 2 $ = 0 (2.1) 

For translatory harmonic excitation in y-direction of the form y(t) = y e*^, 
the boundary conditions are 


9$ 

9y 


= y iOe 
o 


ifit 


( 2 . 2 ) 


at the side walls y = 0 and a 


2 



(2. 3) 



at the container bottom x = -h. The free fluid surface condition is given by the 
kinematic and dynamic conditions, expressing that the normal velocity of a fluid 
particle at the free fluid surface is equal to the normal velocity of the free liquid 
surface, and that the pressure at the free fluid surface is equal to the ambient 
gas pressure, i. e. , p = p . The dynamic condition is obtained from the linear- 
ized unsteady Bernoulli equation and together with the linearized kinematic con- 
dition, yields the free fluid surface condition 


9 2 <i> 

at 2 


+ g 


8 

8x 


= 0 


(2. 4) 


at the free fluid surface x = 0. A transformation 
<My, X, t) = ^0(y, x) + i% o jT|e lfit 


(2. 5) 


yields homogeneous wall boundary conditions. The problem that has to be 
solved now, therefore, is the Laplace equation V 2 0 = 0 for the disturbance 
potential 0 with homogeneous wall boundary conditions. 


30 

9y 


= o 


at y = 0 and a 


80 

ax 


= o 


( 2 . 6 ) 


(2.7) 


at x = -h and the transformed free surface condition 


g 


80 

3x 


fl 2 0 = ifl 2 y y 


( 2 . 8 ) 
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at x = 0. 


A solution satisfying the Laplace equation and the wall boundary condition 
( 2. 6 and 2. 7) is given by 


<M y, x) = Yj 


n = o 


A 

n 


cosh 


J®- (x + h) 

LJL 


cos 




irti ~j 

a J 


cos 



where the values A are unknown coefficients that shall be determined from the 
free fluid surface condition (2. 8). Expanding y into a Fourier cosine series 


y 


a 

2 


OO 


4a y 

- 7]-2 L 


n = 


1 


COS 

2n - 1 

a "L 

(2n - l) 2 


the expression (2. 8) yields 




from which with the square of the circular natural frequency ( see Table 1) 


co 


2 

n 




n = 1, 2, . . . (2. 9) 


(This is obtained for free oscillations of the liquid where y 0 is identically zero. 
The equation (2. 8] is, therefore, satisfied for y Q = 0 and fl 2 = w 2 n(u> 2 ), as can be 
seen from the last equation before equation [2. 9] ) . 

one obtains for the unknown values A^ the expressions: 

4i^ s y a 

A o = ly o^ 2 ; A 2n = ° and A 2n - 1 " 7r 2 (2n - l) 2 (fi 2 - w| n _ ^ 

The disturbance potential is, therefore, 


4 



<M y» x) = -iyjla.( 


1 4 y 

2 " 7T 2 LJ 


00 cosh 


E 


( 2 n - 1 ) 7 T 


(x + h) 


) COS - 


( 2 n - 1 ): 


a 


n = 1 ( O 2 - ( 2 n- i ) 2 cosh ^ 2n ^^ h | 


and the velocity potential yields 



( 2 . 10 ) 


The free fluid surface displacement, the pressure distribution and the velocity 
distribution, as well as the liquid force and moment can be determined from the 
velocity potential by differentiations and integrations. The free liquid surface 
displacement measured from the quiescent fluid surface is given by 


ft 2 m 

x(y * t) =~ y Q e 


( a\ 4a y ^2n- 1 
( y -lj * 2 n V <4 n-i- 1 ) 


cos 


1 


2 n- 1 ) 7 r 


( 2 n - 1 ) ' 


( 2 . 11 )' 


The pressure distribution p = - p34>/8t - pgx is given by 

1)7T 


P = P «V/ !t <(y - f -|r E 


00 cosh 


(2n- l) 7 r , . , , 9 

rT-^ (x+h) J”2„-i 


n -1 ( 2 n- l ) 2 cosh h (I- 17 2 .) 

a ' 2 n-l 


cos 


1 


( 2 n-l)- 


Pgx 
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At the tank walls y = o and y = a the pressure yields 

(2n- 1) 7r 


00 cosh 


(x+h) 


„ 2 liit/a , 4a v 

P y = o = ‘ pS! y o e (2n _ 1)2oosh r M 


) ]’ I 2n-l 


lM h j (1 . 4n _ 


pgx 


and 




lln-l 


( 2 . 12 ) 


n = 1 ( 2n- 1) cos 


J lZn^ih r - 

a 




pgx 


(2. 13) 


At the container bottom x = -h the pressure distribution is given by 


bottom 


= P« 2 y 0 e lflt \( y ■ f ) 


4a y 
- j L 


cos 


E (2n-l)7r I ? 

~r“ y Jv 1 


n= 1 (2n- l) 2 cosh ~ h ^ 1_ ^2n-l|j 


+ Pgh 

(2. 14) 


Integration of the pressure components yields the fluid force and moment. 
The force in y-direction (i. e. , the force per unit width) is 


F = f° (p - p ) dx 

v J , v = a v = o 


which yields with m = pah as the liquid mass per unit width the expression 


6 



F = mfi 2 y Y. 

y ° I n=l <2n- Wf J _ 


00 tanh 


[~!(, 2 n-l)jr 1 2 

L a h _Kn-l 


(2. 15) 


The moment of the liquid about the center of gravity of the undisturbed fluid 
is given by 

M = - f (p -p ) fx + dx - f p (y - dy 

z ' h y = a *y = o V 2/ J Q ^bottom \ y 2/ y 


and yields the expression 


2 i^t 


M = - ma£2 y e 
z J o 


“( 5 ) 


4 2 


00 tanh 


f[2n- IJttJ 

. 2 

2 ± 

l a h ) 

( 7rh/a) (2n- 1) 

_cosh[(2n- l)7rh/a] 


n=* 1 


(2n - i) 3 


V. 


2n- 1 


The velocity distribution is given by 


( 2 . 16 ) 


8# ifit 

v = ~ — = iS2y e 
9y J o 


1+ ±f 4 _ i c ° sh (X4h) [ (211- iw -jl 

% L =i ( i-T in ; 1 )< 2 „-i)co S hp^h] 1 8 L » y J lj 
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which represents the velocity in y-direetion, and 


u = 


9 $ 

9x 


4i£2 

7T 


ifit v 

y~ e L 


4n-l ■toh[^ a2 (x->h)] 


n = i ) (2n-i)cos 


jM h r os 



which is the velocity in x-direction. 

From the results of this section the following is concluded: The natural 
frequency ( 2. 8) of the liquid is inversely proportional to the length of the res- 
ervoir, indicating that these values shall be very small. Since the liquid height 
ratio h/a is very small, the effect of the hyperbolic tangent function is pronounced 
for many of the lower vibration modes, and approaches with increasing mode 
number n the value of unity. The first term in the free fluid surface elevation 
(2. 10) is the surface plane elevation caused by the change of the acceleration 
vector during motion. The infinite series accounts for the liquid waves on the 
surface. In the fluid force expression (2. 15) the first term represents the iner- 
tial force of the liquid. In or close to resonance the liquid force can be a multi- 
ple of the inertial force, therefore endangering the total sound suppressor system. 
In the expression (2. 16) the first term of the liquid moment is identified as the 
moment caused by the shift of the center of gravity for a planar free fluid sur- 
face. In the velocity distribution v the first term represents the container motion. 
Omission of this term yields the velocity distribution in the container. 

As can be seen from the previous results, the liquid motion in the con- 
tainer exhibits singularities at the resonances. However, since there is a 
little damping in the system, (due to internal and wall friction) , a damping 
value should be introduced to limit the amplitudes at the resonances. This will 
be performed with the help of a mechanical analogy, which also will provide a 
simple method for the description of the liquid motion, and reduce the infinite 
number of degrees of freedom of the liquid system to a finite number that is 
pertinent for the dynamic behavior of the liquid. 
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SECTION ill. MECHANICAL ANALOGY OF LIQUID BEHAVIOR 


For the structural analysis a simple mechanical analogon for the de- 
scription of the liquid motion is needed. This mechanical model, by comparison 
with the liquid theory, will then present the amount of sloshing modal masses 
and provide sufficient information of the number of liquid vibration modes that 
must be considered in the structural analysis. In a vibrating container the 
liquid oscillates in the close proximity of the free surface. Part of the liquid 
in the lower part of the container follows the motion of the container like a rigid 
body. To describe this motion with a simple mechanical analogy a spring-mass- 
system is employed. The sloshing mass m^ of the n™ vibration mode is con- 
nected with two springs of stiffness k^^ to the container walls y = o and y = a. 
Its location from the center of gravity of the quiescent liquid is denoted by h 


n 


The non-vibrating liquid in the lower part of the container is represented by a 

This mass is attached rigidly to the container 


mass m of moment of inertia I 
o o 


at a height h^ below the center of gravity of the undisturbed liquid. 


Once the equations of the mechanical analogy are derived, they will be 

compared with the fluid theory to obtain the mechanical values, such as m^, 

m , k , h and h ( Fig. 3) . 
o n n o 


A. Analytical Model for tine Description of the Liquid Motion 


The equations of motion of the mechanical analogy are derived with 
the help of the Lagrange equations. If y^ is the displacement of the n^ 1 sloshing 

mass m^ with respect to the container wall y, the tank displacement in y-direc- 

tion, and i? the rotation about the z-axis, tne kinetic energy is then 

m. 00 

T'l 2 <y + V> + iV >, 4 2 “A + y + V> 2 • < 3 -‘> 

n= 1 


The first two terms represent the kinetic energy of the mass m which is rigidly 

connected with the container. The series describes the kinetic energy of the 

modal masses m . The potential energy is given by 
n 


9 



(3. 2) 


oo oo CO 

V = ^ gt? 2 m h gi? 2 T. mh-gi?y my ky 2 

9. ° o 0 2 z — ' n n ° ' rr n 9 ' rr n 


n= i 


n n 


n = 1 


n n 2 


n= 1 


n n 


Here, the first term represents the potential energy caused by the lifting of the 
mass m during rotation about the z-axis, while the second and third terms 
describe the same fact for the modal masses. The last term is the energy that 
is stored in the springs of stiffness k^. 

To make the mechanical analogy equivalent to the fluid system , the sum 
of the nonsloshing mass and the modal masses must be equal to the total liquid 
mass 


m = m 

o 


OO 




(3. 3) 


Furthermore, for small oscillations, the center of gravity of the liquid shifts 
in the first approximation horizontally only. Therefore, 


m h 
o o 


Y, m h 
u . n i 
n = 1 


(3. 4) 


must be satisfied. Because of this condition, the first two terms of the potential 
energy cancel each other. 

The equations of motion are now derived from the Lagrange equation 


_d 

dt 





(3.5) 


where Q. are the generalized forces and are 


Q = -F ; Q = -M ; Q - 0 

y y * 2 y n 


(3. 6) 
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and y, ■& and y are the generalized coordinates. The equations of motion are, 
therefore. 


mj'j + hj) + Yj m n (y n + y +h n 2) =-F 

n= 1 


(3. 7) 


I ^ + m h 2 3 - g Y m h + Y m „ h n ( yn +h n^ ) = - M „ ( 3. 8) 


o oo “ L -‘ . nn t - J , nnn n 
n=l n=i 


y + w 2 y = - y - hi? + gd- ; n=i, 2,... (3.9) 

nnn n 

k 

where w 2 = (see 2. 8). The first equation is the force equation, the second 

n m 

n 

is the moment equation and the third represents the sloshing equations (n = 1, 

2 , . . . ) . 

B. Solution of the Equations of the Model and Determination of the 
Mechanical Analogy Values 

Since a pitching motion (i.e. , oscillation about the z-axis) of the 
sound suppressor reservoir is unlikely, only the translatory excitation case is 
treated. If S- = 0, the equations of motion yield 


my 


y + Z 


n= 1 


n n y 


(3. 10) 


Y m(hy - gy)=-M 

n n n n z 

n= 1 


(3.11) 


y* n + «^y n = - y ; (n = l, 2, . . . ) 


(3. 12) 
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iQt , _ lfit 

From the sloshing equation with y = y^e and y^ = Y^e 




y — — — 

n co 2 - S2 2 
n 


n 2 y 

° r y n = J^ 


n 


(3. 13) 


is obtained. Therefore, 


»• I2 2 y 

y n _ co 2 -0 2 . (3.14) 

n 


With this, the force equation yields 



The moment is then 


M z = -my 2 
n = i 


m (g + h S2 2 ) 
n & n 

m (cj 2 - S2 2 ) 
n 


(3. 15) 


(3. 16) 


After this derivation and solution of the mechanical analogy, the mechan- 
ical values have to be determined. This will be performed by comparing the 
mechanical analogy with the results of the liquid theory (Section II). However, 
the moment equation (2. 16) as obtained from the fluid theory is not compatible 
with the moment equation of the mechanical analogy (3. 16) and has to be trans- 
formed. 


The other expressions are then compared. First, the spring- stiffness 

th 

k has to be chosen such that its ratio with the n sloshing mass m represents 
n n 

the square of the circular natural frequency: 


k 

_J1 = ^2 = SIZE 

m n a 
n 



for n = i , 2, • , . 
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Comparison of the forces (2. 15) and (3. 15) yields (Table 2) 


, , (2n-l)7r, 

21-1 a tanh [_-~r ~ h J 

m 7r 3 ( 2n- 1) 3 (h/a) 


for n - 1, 2, . . . 


(3. 17) 


and 


m 

With 

oo 

1 A V — 1 

12 (h/a) 7r 4 ^ ^ (2n-l) 4 (h/a) 


the transformed moment reads 


= - my 


V 8a 


taih [ (2p a ~ 1 > ,rh ] 


n ^ 1 7T 4 (h/a)'(2n-l) 4 (w| n _ 1 -fi 2 ) 


g( 2n-l) 7 r 
a 


+ fl 2 \~ (2n-l) - 2 tanh 


(2n-l) 7r 


2a 


(3. 13) 


Comparison with the moment equation (2. 16) and that of the mechanical 
analogy yields the relation for h 

n 


h = h (a _ 4a 
n 2 7rh(2n-l) 



(3. 19) 


The magnitude of the modal masses and their location depends only on 
the geometry of the container. Therefore, it will be expected that the influence 
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of the liquid in the reservoir upon the structural behavior of the sound suppres- 
sor system can be determined mainly by the geometry of the reservoir. 


SECTION IV. STRUCTURAL ANALYSIS 


For the proper design of the system, the knowledge of the interaction of 
the structure of the sound suppressor system and the liquid in the reservoir is 
of great importance. This is particularly true since more than eighty percent 
(four million gallons) of the liquid sloshes in the container. In the analysis the 
liquid in the container is described by the mechanical analogy ( see Section III) . 
Since the occurrence of longitudinal (in x-direction), pitching, yawing and roll 
motions of the container is very remote, analysis is restricted to the tranlation- 
al degrees of freedom in y and z direction. These, however, can be assumed as 
uncoupled. Thus , a two degree of freedom system is obtained if all higher 
sloshing masses are neglected and only the fundamental sloshing mode is re- 
tained. This is justified since the second sloshing mass is only 8. 6 percent of 
the total liquid mass. Since the mass of the support structure exhibits a rather 
large portion of the total mass of the system, it cannot be neglected. The amount 
of 17/35 of the support structure is added to the mass of the container structure 
m and the nonsloshing mass m because the support structure represents 
cantilever beams. The container is carried by v pairs of cantilever beams on 
either of the long sides in the shown orientation ( Fig. 4) and by p pairs of 
cantilever beams at the small sides of the container (z-direction) . The stiff- 
ness of these beams are El in y-direction and El in z-direction. The beams 

y z 

are supposed to be identical except for the orientation, which is exhibited in 
Figure 4. The value E is Young's modulus of elasticity, I is the geometric 
moment of inertia of the cross section and I is the length of the beams. The 
sloshing mass of the n^ 1 liquid mode is attached with springs of stiffness k /2 

and k /2 to the container walls at a height h . The distance from the grade is 
zn n 

denoted by L^. The mass M composed of the container mass m^, nonsloshing 

mass m and 17/35 of the support structure, i.e. , 
o 

M = m + m + 17/35m 

co ss 

is located at a distance l from the grade. Its displacement in y-direction is 
denoted by y^ and in z-direction by z^. The n^ slosh mass performs a 


14 



displacement relative to the reservoir wall and is denoted by f° r motion in 

the y-direction and by z^ 2 for motion in the z-direction. If, as in the numerical 

evaluation, only the fundamental slosh mass is considered, its displacement 
relative to the container wall is denoted by y 2 for the y-direction and by z 2 for 
the z-direction. 

A. Equations of Motion 

Since the fluid mass that participates in the motion of the free fluid 

surface (i. e. , the sloshing mass) decreases as l/(2n-l) 3 , only one sloshing 

mass m = m, attached at a location L = I + h + h . from the lower end of the 
1 o 1 

beam, is considered (see Fig. 4). The sum of the mass of the container struc- 
ture, the mass of the non- sloshing liquid and the effective mass of the support 
structure is denoted by M. The generalized coordinates are y . and y , where 

y^ is the displacement of the mass M at the location i and y 2 is the displacement 

of the sloshing mass m at the location L. For the derivation of the equations of 
motion the displacement influence coefficient method is employed. For this 
reason the displacement of the beam at the location x caused by a force F at the 
location £ has to be determined. The moment at the location x caused by a force 
at the location £ is given by 


M = Ely " = F <£ - x) for x < £ . 


In assuming a uniform beam (i. e. , constant stiffness El) , the displace- 
ment can easily be obtained by integration of 


y"=ti <*- x) 


for x < £ 


with the boundary conditions of zero slope and zero displacement at x = o in the 
case x < £. This yields a slope (x < £) 



and a displacement for x < £ 
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y (x) . 


For x > £ the displacement is 
y(x) = y(£) + y’(|)(x- |) 
which yields (x > £) 

y(x) = ^{^ x2_x3+ (x ‘ ' 

The influence coefficients S„ are determined next. Applying a unit force 
at the location ( 1) yields at this location with x = | = £ a deflection 


6 ^ - 


I 3 
3 El 


(4.1) 


The deflection of the mass m (i. e. , at the location [2]) is then with x = L and 

« =i 

s «=iilf 3iL 


which yields 


- L s + ( L - i ) 


■ 



A unit force applied at the location ( 2) , i. e. , at the mass point m, yields 
a deflection of that point by 

_ kL 3 + 3EI 

° 22 3kEI (4.3) 
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I 


This is obtained from the fact that the springs with stiffness k and 3EI/i 3 are in 
series, which means that their influence numbers add. The equations of motion 
for the total system are now derived. 

A unit force applied in y-direction at the location ( 1) , i. e. , the location 
of the mass M, yields at that location (1) a deflection 6 ^ ^ in y-direction and 

a deflection S L in z-direction. At the location ( 2 ) , i. e. , the location of the 
Vl z i 

sloshing mass, the displacement is 6 ^^ in y-direction and 6 ^ ^ in z-direction. 

A unit force applied in z-direction at the location ( 1) yields the deflections 

<5j in y-direction and 6 , in z-direction at the location ( 1 ) , and 6-, in 
ZlYl z i z i lz iy 2 

y-direction and in z-direction at the location ( 2 ). Similar results for the 

z l z 2 

deflections are obtained by applying a unit force at the location ( 2 ) first in y 

and then in z-direction. Therefore, &> is the deflection in y-direction at the 

y 2 yi 

location ( 1 ) caused by a unit force in y-direction, applied at the location ( 2 ). 

62 is the deflection in z-direction at the location ( 1 ) caused by that unit force, 

y 2 z i 

<5 2 is the displacement in y-direction at the location ( 2 ) and 6 ? is the de- 

y 2 y 2 y 2 z 2 

flection in z-direction at the location ( 2 ) caused by a unit force in y-direction at 
the location (2). A unit force in z-direction at the location (2) yields a deflec- 
tion (5 2 in y-direction and 6 ? in z-direction at the location ( 1 ) and 6? 

%yt Szi z 2 y 2 

in y-direction and 6 2 in z-direction at the location ( 2) . The static problem 

z 2 z 2 

is therefore solved if the influence coefficients 6. and 6 . (i, j, k = 1 , 2 ) 

are known. The deflections caused by the forces F 1 = Y 1 j -t- Z j k at the location 
(i) and F 2 — Y 2 j + Z 2 k at the location (2) are, therefore: 


yi - « 5 i Yj + 6j Z i + d 2 Y 2 + 60 Z 2 

Yiyi z iyi y 2 yi z^i 


Zi — 6 -i Y< + 6j Zi + d 2 Yo + 6o Z 2 

Yi z i z i z i y 2 Zi z 2 z t 
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Y2 = 
z 2 = 


61 Y t + 6i 

W 2 z iy2 

6 t Yi + 6i 

Vl z 2 z l z 2 


Zi + 6 2 
Zj + 62 


Y2Y2 


Y2 z 2 


Y 2 + o 2 
Y 2 + 62 


z 2Y2 


^2 Z 2 


z 

z 


2 


2 


Since the beams are all parallel they act like parallel springs, thus having a 
spring stiffness of the sum of their stiffnesses. 


Considered is one pair of beams of the system. If its geometric moment 
of inertia for bending in y-direction is given by I and in z-direction by 1^, the 

stiffness in y-direction of 2 p pairs of beams on the short sides of the system 
and 2 v pairs of beams on the long sides of the container is given by: 


3 EI 3 El 

j)3^ 2p 


The stiffness in z-direction yields: 


3 El 

Y 

£ 3 


3 EI 

2 m + -jr- 


2v . 


For free oscillations the forces acting on the system are 


Yj = - My\ , Z t = -Mz\ 

Y 2 = -my 2 , Z 2 = -mz 2 


and the equations of motion are given by 


Yl + 

z l + 

Y2 + 
z 2 + 


\ T _ My*i + 61 Mz'j +62 my 2 + o 2 mz 2 = 0 
TlYl z lYl Y2Y1 z 2Yl 


61 My. + 6i Mz*i + 60 my, + 6, mz, = 0 

Vl z l 1 lz l z l 2 y 2 z i 2 2z 2 z i 


My\ + 61 Mz'i + S 2 my 2 + S 2 mz 2 = 0 
W2 Z U2 YtYl Z 2Z2 


6 Myj + 61 Mz'i + < 5 2 my 2 + o 2 mz, = 0 
yi z 2 1 z x z 2 y2 z 2 Z 2 Z 2 


( 4 . 4 ) 


( 4 . 5 ) 


( 4 . 6) 


( 4 . 7 ) 
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Since the support beams are chosen to be I-beams , thus exhibiting sym- 
metry, the influence coefficients 61 , <5i , S 2 > . <5i , <5 2 , 

* z iyi ti z i z 2yi T2 z i lz iy2 ^2 

61 , and 60 vanish. This means that a force in y-direction does not create 

j 1 z 2 j 2 z 2 

a deflection in z-direction. Mechanically this means that the coupling of the y 
and z-direction is zero. The equations of motion, therefore, split up into two 
pairs of equations and are 


y, + + Vi ” 5 ’ 2 = 0 




(4.8) 


and 


Zi + <5i Mz'i + 60 mz*o = 0 
1 1 z 1 z I i z 2 z i 


Zi + 61 Mz'i + 60 mz 2 = 0 
zjz 2 ‘ J z 2 z 2 


(4. 9) 


The influence coefficients with (4. i) through (4. 3) are given by 

i. <sl 

. 


Viyi 6E(vI +jul ) 

y z 


i' 


#■] 


Viy2 12E(i'I +/il z ) 


,11 O) cti 


, <52- 


V 2 yi i2E(vi + fxi ) 
y z 

k L 3 + 6E( vl + \x I ) 

= _y y z 

y2y2 6kyE(pI y + Ml z ) 


(4. 10) 


and 


z i z l = 


6E( v\ +ixl ) 

z y 


i® 


[# ■] 


i z i z 2 l2E(i^I + ix I ) 
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6 


h 2 Zi 



12E( vl + 
z 


'‘V 


k L 3 + 6E(vI + Ml ) 

6 z__Z_ 

2 z 2 z 2 6k E(vl + mI ) 
z z y 


(4.11) 


As can be seen from the previous equations, the two pairs of equations of 
motion can be solved independently. Transformation of the displacement influ- 
ence coefficient matrix 


[6] = 


’V 


u 2 

lYl T2Y1 


6l yiy2 62 y2y2 


i • 


6vEI 


# *] 


12vEI 


*(t)- 


12vEI 


k L 3 + 6vEI 

_Y Y 

6k El r 

y y 


into the force influence coefficient matrix yields with k^= (-1) 


A 

A + v v\ 


where 


A is the determinant of the displacement influence coefficients and A the 
subdeterminant to the term (\v ) . 


[K]=J 


24EI p{k L 3 + 6EI v) 

y y y 


■12I/EI k 1 3 

y y 


3 7-1 


■12i/EI k f 3 

y y 


(t) 


241 3 vEI k 

y y 


kil k 12 
k 21 k 22 


J 


It is 


A = 41 3 (k L 3 + 6^EI ) - l 6 k 3 fy) - 1 

y y y L \* J 
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The equations of motion in y-direction, therefore, can also be written as 


My t + k^y! + = 0 


my 2 + k 21 y t + k 22 y 2 = 0 . 


These equations exhibit static coupling. 


B. Analysis of the Interaction of Structure and Liquid in 
the Reservoir 


Since the equations of motion in y and z-direction are the same kind, 

the solution of one set is sufficient. The final results can be modified by writing 

just the value z (i = i, 2) instead of y., and the expressions 6 . instead of the 

i i iy k 


influence coefficients 8 . 


Vk 


1. Free Oscillations. 

For the free oscillations the solution is 


Yl = 


y2 = 



cos cot 
sin cot 

coscot 
sin cot 


} 

) 


which yield the frequency equation 


1 - 6 l Mco 2 

Viyi 


1 - 6i Mco 2 
W2 


- 8a mco 2 
Y2Yl 


1 - 6 ? mco 2 
Y2Y2 


= 0 
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This yields with 1/co 2 = X 2 the expression 


X 4 - X 2 


{%iyi 


M + d 2 m 1 + mM)6i d 2 - <5i 6 2 

y 2 Y 2 J \ TlYl Y 2 Y 2 TiY 2 Y 2 Y 1 




The solution of this equation is 


X 2 =±(d t M + S 2 m\± 7 : 
1/2 2 ^VlYl T2Y2 




M - m6 2 ) + 4mM<5- L 6 2 ) 

y 2 Y2 TlY2 YzYl 

(4. 12) 


which gives the square of the natural circular frequencies 


on 


i = i/Xi 2 and co 2 2 = 1 /X 2 2 • 


The Amplitude ratios are 


An 6. 


11 °2 va mwi 

_ ypy?. 


mui 2 - 1 


An 




A 12 62, 


_ y 2 X 2 

* 22 


mw 2 2 - 1 1 


A 99 d-i Mw 2 2 


X 2 


(4. 13) 


These ratios are equal to the ratio of the displacements yi/y 2 at these frequen- 
cies. With these functions 


yn (t) - A^ cos co jt 


, yi2 (t) = A 12 cos co 2 t 


y 21 (t) = AjA^ cos coit 


y 2 2 (t) — ^A^ cos OJ 2 t 


are obtained, which are particular solutions with the integration constants A tl 
and Aj 2 . A solution of the form 
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yj = Bj sin o>t 
y 2 = B 2 sin cot 

yields the same frequency equation and mode shape numbers Aj and Ag. The 
complete solution, therefore, is given by 


yj(t) = A tl cos cojt + Bjj sin cojt + cos w 2 t + B 12 sin w 2 t 


v»( t) = A.f~A 


where Ajj, A^, Bn and B^ are four integration constants which can be deter- 
mined by the initial conditions at t = o. With 


Yi(o) =y 10 , yi(o) = u 10 

y2<o) =Y20 » y 2 (o) = TJ20 (4.14) 

the solution is 


yi(t) 



(4. 15) 


A i f 1 1 

y^t) = (yi(A 2 - y 2 o) cos Wit +— ( vnjAz - D 2 o) sin Wit 

** r i i 

+ ( A 2 - A 4 ) * y20 " Ajyi °) COS W2t ~co 2 * Vzo " Xl ^ Sin Wzt 


(4.16) 


Similar results are obtained for the free oscillation in z-direction. 
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2. Forced Harmonic Excitation. 


For a forced oscillation with an excitation function F sin fit 

_ oy 

z 

acting at the container, i. e. , at the location ( 1 ) , the equations of motion are: 


y ‘ + V, My > + V myi * V.V sin nt 


y! + V“ y ‘ + V ” y2 “ V. r oy si ” nt 


(4. 17) 


for forcing in y-direction, and for a forcing function in z-direction at the loca- 
tion ( 1 ) it is 


z, + 64 Mz, + do mz 2 = 64 F sin fit 
Z]Z X z 2 z l z l z l ° z 


z ? + di Mz, + & mzo = 64 F sin fit 

ZjZ 2 1 Z 2 Z 2 z l z 2 OZ 


(4. 18) 


The steady state solution of the system (4. 17) yields the response amplitudes 


<5, [ 1 - mfi 2 <5 2 J + <$4 d 2 mfi 2 

Y = v 34Y1 Y2Y2 y.iy?. y?.yi 

1 oy mMd, d 2 (fi 2 - u t 2 ) (fi 2 - u> 2 2 ) 

Tiyi y2y2 


(4. 19) 


Y, = F 


m 


oy mMd^^d^ _ (fi 2 - o^ 2 ) (fi 2 - co 2 2 ) 


Y2Y2 


(4.20) 


Similar results are obtained for the system (4. 18) by just exchanging y. and Y. 

with z. and Z . 11 

1 1 
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3. Wind Response . The response of the system due to some wind 
loading is important for the proper design of the structure. For this reason, 
the effect of an arbitrary forcing function F (t) and F (t) in y and z-direction 
respectively is investigated. ^ 


The equations of motion for excitation in y-direction are given by 




yz + 6 V 1 y ! My ‘ + 6 w*« ' v, v * 1 


and for excitation in z-direction, the equations of motion are 


Zh + Mz'i + (5 2 mz 2 = 61 F (t) 

1 1 ZiZi z 2 z i lz l z i z 


Zo + 61 Mz’i + 6n mz 2 = 61 F (t) 

lz iZ 2 ^ Z 2 Z 2 z l z 2 z 


(4.21) 
(4. 22) 


(4. 23) 


(4. 24) 


These represent two linear inhomogeneous differential equations, of 
which the transient solution shall be determined for an arbitrary excitation. 
Since the system already might be in some kind of motion caused by the surge 
wave or some disturbance created by the firing of the rocket engines before the 
wind hits the system at a time t = t Q , the free oscillations during the time inter- 
val 0 < t ^ t Q are also treated in order to present properly the initial conditions 

at the time t . 

o 


a. Step function. For the derivation of an arbitrary time function 
input, the case of a suddenly applied constant force has to be treated first. Since 
the above systems of two differential equations are of the same form except for 
the constant coefficients, the derivation of the response is restricted to one sys- 
tem. Introducing for 

6-l = < 5 , , do =61 = & 2 and 6 2 = d 2 

Tiyi 1 y2Yi VU2 Y2Y2 
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the equations of motion that have to be treated for a forcing function ( Fig. 5) 


F(t) =< 


0 for 0 ^ t ^ t 


F for t 2 : t 
o o 


are given by 

yi + SjMyj + fi^mya = <5iF(t) 
y 2 + S^My'j + 6 2 my 2 = 6 12 F(t) 


(4.25) 


(4. 26) 
(4. 27) 


For the time interval 0 < t t , the system is not disturbed and exhibits zeros 
on the righthand side. With the initial conditions at the time t = 0 


Yl(o) = yio 
y 2 (o) = y 20 


and 


yi(o) = v 10 
y 2 (o) = v 20 


the solution yields for the time interval 0 ^ t ^ t Q : 


yi = 


( A 2 


^rf 


(yi 0^2 “ Y 2 o) cos + ( v 10^2 - V 20 ) si 11 


+ (y 2 0 - ^-lyio) cos w 2* + 7 T ( v 20 - X l v io) sin W 2 t 

OJ2 


•) 


y2 = 


(a 2 


- M) { ! 


Ai(yioX 2 - y 20 ) cos o^t + (v^ - v 20 ) sin cojt 


wi 


(4.28) 


+ A 2 (y 20 - Atfio) cos co 2 t + “*■ (v 20 - XiV 10 )sin 

a>2 


in co t L 

2 J ‘ 


(4.29) 
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At the time t = t Q , a constant force F^( wind) is suddenly applied. The solution 

of the above differential equations in this time interval (t = t ) is 

o 

yi = cos GOit + B tl sin cojt + A^ cos o> 2 t + B 12 sin <u 2 t + ^ 

t > t 
> 0 

y 2 = Aj( A tl cos wj t + B 1A sin u>*t) + A 2 (A 12 cos o> 2 t + B 12 sin w 2 t) r 

+ d^F > 

u o 


where A^, A 12 , Bjj and B 12 are integration constants and ^F^ and 6^1^ are the 

particular solutions connected with the inhomogeneous part of the differential 
equations. The constants of integration can be determined by the new initial 
conditions at t = as obtained from the final values of displacement and velocity 

in the first time interval. This means that at the time the wind disturbance ex- 
cites the system, it had already some displacement and velocity from a previous 
disturbance. The constants of integration finally are 


A U 

B« 

A 12 
B 12 


(yi0^2 - Y2 o) 

F Q (^1 A 2 - ^ 12 ) 

(A 2 - Aj) 

(A 2 - Aj) 

v 10^2 “ v 20 

F q ^1 A 2 ~ ^12) 

w l(A 2 - Aj) 

(A 2 - Aj) 

Y20 ~ Ajy 10 

F Q ( 6 12 “ A l fi l ) 

(A 2 - Aj) 

(A 2 - Aj) 

V 20 - *iV 10 

F Q ( 6 12 “ Ajdj) 

~ Aj) 

(A 2 - Aj) 


(cos WlO 


(sin w it ) 


(cos w 2 t Q ) 


(sin o^t^) 


and the solution of the system yields for 


t > t 
o 
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(y 10^-2 ~ 72 o) v 10^2 - v 20 Y20 “ Xjy 10 

y< = 7 — cos wit + 77 — r sin w it + — — cos w 2 t 

J1 X 2 - Xj 1 wj(X 2 - Xj) X2 - Xj 


v 20 " Xtv 10 
w 2 (X^ - Xj) 


sin co 2 t + djF^ 


~ ^12) 
X 2 - Xj 


COS co j( t - t Q ) 


F ( 6^2 ~ Xj6j) 
X 2 - Xj 


cos w 2 (t - 


V 


(4. 30) 


Xi(yioX2 - y2o) 
y2 ” x 2 - Xj 


Xi( v 10 X 2 - v 20 ) 

cos Wit + — 77 r-r — sin o^t 

1 - X t ) 


X 2 (y 2 o - XiYio) 


X 2 (v 2 o - XjVjq) 


X 2 


— cos + ^ z Xl) sln + 6 «= F o 


fXiF (6iX 2 - 6^) ^ C X 2 F (5j2 - 

-l - — 008 “ l<t - ‘o !/ - 1 


- X^i) 


cos co 2 (t 

(4.31) 


-*o} 


It can be seen immediately that if the system was at rest at t = 0, i. e. , y 10 = 
y 2 o = 0; vio = v 2 o = 0, the motion for t s t Q is given by expressions in the second 

line of each formula. It represents oscillations with the circular frequency w* 
and w 2 about the new equilibrium positions SjF^ for y* and S^F^ for y 2 , indicating 

that the structure as well as the free fluid surface oscillate about an inclined 
position. 

If the system is at rest until at a time t = 0 a wind step function of mag- 
nitude F is applied to the container, the system responds with 
o 

( <6iX 2 - <5^) (<5^ - Xjdj) 

yi = F S — — cos wjt - — cos w 2 t 

° L X 2 - Xj X 2 - Xj 




( 4. 32) 
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Yt 


- F „ - 


^l( - 612) 

^2 " 


COS OJjt - 


^2(^12 - ^-i^i) 


cos w 2 1 


}. 


(4. 33) 


b. Rectangular pulse. With these results a rectangular pulse of 
magnitude F q can be obtained by subtracting the response of a step function ap- 
plied at t = tj and of magnitude F . It is, therefore, the response of the system 
due to a rectangular pulse of magnitude F q and duration tj seconds given by the 
equation ( 4. 32) and ( 4. 33) for 0 ^ t ^ tj and for t 2: tj by 


Yi = 


F 

o 


{ 


<5 1^2 - <5i2 
(A 2 - Ai) 


[cos CO t - tj) - COS CO it] 


( S 12 - Ai6i) 

(A 2 - Ai) 


[cos co 2 (t - ti) - cos co 2 t] 


■} 


Y2 = 


F 

o 


{ 


Xi( djA-2 - <5^) 
(A 2 - A.j) 


[cos u>i (t - tj) - cos Wit] 


(4. 34) 


A 2 ( dj2 - A-idj) 
A 2 - Ai 


[cos w 2 (t - t^ - cos w 2 t] 


) 


(4. 35) 


Several elementary forcing functions shall be treated. First, the response 

— Q/ 1 

of the system to a simple exponential asymptotic step F ( 1 - e ) shall be 

° -St 

investigated. Then an exponentially decaying pulse of the form F^e will be 
treated, where the latter one represents for appropriately large F q very closely 
the blast loading case. 

c. The exponential asymptotic step. The equations of motion for 
the system disturbed by an exponential asymptotic step of the form (Fig. 6) 
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F(t) F ( i - e ) 

o 


; for ostst 


; for t s t 


(4. 36) 


are given by 


Yl + SiMy*! + b^mfo = \ 6 jF o ( 1 - e “ * 1 V ) 


; for o ^ t ^ t 


; for t s t 


o J(4. 37) 


; for o < t ^ t 


y 2 + S^My'i + 6 2 m y ‘2 = \ 6 12 F 0 ( 1 " e ^ ^ ) 


; for t > t 


o - y (4. 38) 


The solution in the interval 0 ^ t ^ t is given by (4. 28) and (4. 29) if the initial 
conditions at t = 0 are given by yj( 0 ) = y 10 , y 2 ( 0 ) = y 20 , yj( 0 ) = v i0 and y 2 ( 0 ) = 
V 20 . At the time t = t the exponential asymptotic step commences. For t ^ t Q 
the solution, therefore, is 


Yi = An cos w^t + sin cujt + Aj 2 cos cu 2 t + Bj 2 sin u> 2 t + bjF^ 

F 0 {^i 2 2 ma 2 “ <5i(l + <5 2 ma; 2 )^ e 

+ ^1 + SjMo! 2 + 6 2 mo! 2 + mMa 4 (6i6 2 - 6 j 2 2 ) y 

y 2 = A^Ajj cos cojt + B^ sin cojt) + A 2 (A 12 cos a^t + Bj 2 cos co 2 t) 

r- -a ( t - t ) 

I 6 O 

+ fi « F 0 ^1 + <5jMQ! 2 + (5 2 m« 2 + mMa; 4 (<5i<5 2 - b^ 2 ) 
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where A 11( B 11# A 12 and B 12 are integration constants and the last two terms in 
each equation represent the particular solution of the system. The constants of 
integration can be determined by the new initial conditions at the time t = t^. 

These initial conditions are the final values for the displacement and the velocity 
in the first time interval. The constants of integration finally are 

y2o - *iyio cos w 2t 0 
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and the solution of the system yields for t > t the expressions: 
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a^mo 2 ~ 6j(l + <5 2 ma 2 ) 
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(4.39) 
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| 6l2 2 mG; 2 - 6^1 + ^mce 2 ) C 
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( 4. 40) 


For a = 0 the differential equations exhibit no excitation and the solution is that 
of the free oscillation. For small values a the forcing function approaches the 
value F q very slowly, while with an increase in magnitude of the excitation func- 
tion it approaches more and more the step function. For a — °° the previous 
results should, therefore, exhibit in the limit the results of the excitation of the 
system with a step function. 

It can be seen immediately that if the system was at rest at t = 0, the 

motion for t ^ t is given by the expressions (4. 39) and (4. 40) if the first line 
o 

is omitted. 

d. Exponentially decaying pulse. The equations of motion for the 
system disturbed by a wind which is described by an exponentially decaying pulse 
starting at the time t = t with an amplitude F q ( see Fig. 6) are given by 
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if 0<tst 


y 2 + 6i2 M y’l + 6 2 my 2 = 


- R(t - t ) 
6 12 F e PK o’ 
1 o 


if t > t 


(4. 42) 


The solution in the time interval 0 t ^ t is given by (4. 28) and (4. 29) if the 
initial conditions at t = 0 are given by yi(0) = y 10 , y 2 (0) = y 20 , yi(0) = v 10 and 
y 2 (0) = v 20 . The solution of the system yields for t 2 the expressions 
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+ 5 12 
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( 4. 43) 
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For 0 = 0 this yields the response of the system due to a step function at 
t = t^. For increasing 0 the excitation function approaches the value zero more 

rapidly. If the system was at rest at t = 0 the motion is described by the above 
expressions by omitting the first line. 

e. Arbitrary wind build-up. The results of these investigations 
can now be used to describe an arbitrary, suddenly applied case of excitation. 

If at a time t = t a.force F(t) is suddenly applied, the response of the system 

can be found by the superposition of all effects due to an infinitesimal change in 
the force (cross-hatched area in Fig; 7) . The response of the system due to 
tills part is the same as that of a step function suddenly applied at t = t and of the 

dF 

magnitude ~ dr = F'(T)dT. By integration, the effect of the total change in F(t) 

from t = t on is obtained. In addition, the constant part due to F(t ) must be 
o o 

superimposed. The problem that has to be solved now is given by the system of 

differential equations 
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(4. 45) 


(4. 46) 


The solution of this system in the time interval 0 ^ t ^ t is the same as 
that obtained previously ( 4. 28) and ( 4. 29) . For t > t the solution is composed 
of the various effects of the strips F*(r)dr in addition to that of the constant force 
F(t Q ). The response due to F’(r)dr is given by 
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Superimposing the response of the step function of magnitude F(t^) with 
the integrated parts of the force F(t) above F(t^) yields the total response 
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and 


Y2 = 


My 10*2 “ Y 2 o) 
(*2 - 


COS CO it + 


M VioA 2 - V 2 o) 
COi(X 2 “ M 


sin co it 


^(y20 “ Mio) M V 20 - ^i v io) 

+ — cos " 2 * + sin " 2t 


+ F(t ) < 6* - 


Ai( 6 iA 2 - 612) 
(Aa - \i) 


COS a>i(t - t ) - 


A 2 (6i2 - djAi) 


/ /-v -k \ cos w 2 (t-t )> 

o (A,, - A t ) o( 


39 



A-i(5l^2 “ 5 j 2 ) 

(\ 2 “ ^-l) 


COS COi(t - t) 


^ 2(612 " ^l^i) 
U 2 - At) 


- t) } 


cos co 2 (t - t) ) F'(t)c 1 t 


Integrating by parts yields then 
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(4. 50) 


If the system was at rest at the time t = 0 the response for t s t is given 
by the above formulae by omitting the first four terms in each equation. 

SECTION V. DAMPED VIBRATIONS OF THE SYSTEM 


Especially in the harmonically forced vibration case the displacements of 
the sound suppressor can be rather large around the natural frequencies of the 
system. Since the undamped case exhibits singularities at the natural frequencies, 
damping must be included to find the maximum amplitudes. 

A. Introduction of Damping 

The results of the previously derived potential theory of fluid oscil- 
lations are only valid and applicable if the exciting frequency is not too close to 
the resonances and if the excitation amplitude is not too large. The latter con- 
dition can be assumed to be satisfied; thus, linearized theory is justified. In 
the vicinity of the resonances, however, fluid forces occur which are a multiple 
of the inertial force of the liquid. These areas, of course, represent the impor- 
tant frequency range in which the motion of the liquid will have its most pronounced 
effect upon the structure of the sound suppressor system. The- introduction of a 
damping factor is , therefore , of importance for the determination of the liquid 
forces in these frequency ranges. It is introduced as a linear viscous damping 
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of which the damping factor is determined by experiments. This actually repre- 
sents an equivalent linear damping, which, in the case of the sound suppressor 
system, is of very small magnitude ( £ = 0. 001 to 0. 01). 

To account for damping in the liquid a dissipation function 


oo 



c y 
n 


2 

n 


has to be introduced in the description of the liquid motion by the analytical model. 
This yields an additional term in the sloshing equations (3. 9) which with 


c = 2£ cj m 
n n n n 


now reads 


2 „ - 


y +2£u;y + cr “y 
J n n n n n n 


-y 


(5. i) 


and yields a liquid force 


F = my 
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oo 


n = 1 


m 
n 
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g 2 

(c o 2 - £2 2 -f- 2iJ co fi) 
n n n 


(5.2) 


and a moment of the fluid 


M z = ~ m y Z 

n= i 


m 
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(« + V‘) 

(cu 2 - + 2if co £2) 

n n n 


(5. 3) 


After introduction of the linear damping terms into the mechanical model, the 
values for the damped liquid force and moment can be obtained by merely intro- 
ducing in the resonance terms instead of (co 2 - ^ 2 ) the expression (co 2 - £2 2 + 

n n 

2i£ w fi). The structural damping can be introduced as a linear viscous damping 

and is represented as a dashpot with a damping coefficient c at the location of 

s s 

the mass M. 
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Therefore, the equations of motion (4. 8) are with damping, in y-direc- 


tion 


(y) 


(y) 
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ss 
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ss y 2 y 2 


SS 


r t 
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(5. 5) 


and similarly, in z-direction 
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r** . ' • 


(z). 


z 1 + (Mzj + c zj) + <5 2 (Mz 2 + c z 2 ) = 0 (= <5j F sin fit) 


z 1 z 1 ss Z 2 Zj 4 ss 
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(5. 6) 
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z 2 + <5i (Mzj + c zj) + <5 2 (Mz 2 + c z 2 ) =0 (= 6i F sin fit) 

° ±r 7 J r 7 ~ 1 oc 1 c r/~rr~ L o o L *-r, .r, ^ r\rr 1 


zjz 2 ss A z 2 z 2 ss 


z t z 2 oz 


(5.7) 


For the free oscillation problem the right hand sides of these equations are zero, 
while for the forced vibration problem ( excitation again acting at the container) 
the right hand sides of these differential equations are represented by the ex- 
pressions in the parentheses. 

The damping coefficient c of the structural damping is given by 

s s 

c = g Mu = Mg VISt 
ss ss ss ss V 


which in y-direction is 
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c ss s ss M V m 3 


(5. 8) 
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and in z-direction is 


(z) 

C ss = S ss M ’ M^ 3 


(5. 9) 


The value g is the structural damping and probably has a value in the range 
ss 

0. 005 < g < 0. 02. 
ss 

The damping coefficient c of the sloshing liquid is given by 

s 


c = 2mw £ 
s s s 


which in y-direction is 


c = 2mco £ 
s sy s 


and in z-direction is 


c = 2mw £ 
s sz s 


(5. 10) 


(5.11) 


B. Damped Sinusoidal Response 

To obtain the steady-state solution for the damped system caused by 
some harmonic excitation, complex solution methods are employed. With 

di = 6a, do = <5.io and 6o = 6o, the two simultaneous differential equations 
Viyi yjy2 T2y2 

( 5. 4) and ( 5. 5) can be written in the form 
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(5. 12) 
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(5. 13) 
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The steady-state solution of this system is 


yi = Yie 
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Y 2 = Y 2 e 
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which inserted in the above equations yields two complex algebraic equations for 
the complex amplitudes Yj and Y 2 : 
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Since the amplitude of the structure and liquid rather than the phase relation is 
of interest, only the absolute value of the amplitudes is investigated: 
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A parameter a has been included to show the effect of stiffness variation. 
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C. Damped Liquid Oscillations 


The introduction of linear damping into the results of the liquid anal- 
ysis (Section II) requires in the resonance terms the value (i - ?7 2 + 2i £ rj ) 

2 

instead of (1 - tj )* Some of the expressions in Section II, however, cannot be 
used in their present form. A transformation with 
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( 2n-l ) 

. „ 2 m 00 cos ttx 

4aft^ iS2t v La 

x(y t) = — 5 — Y e /, 

J » n l s o u . (v 




« ° n = i ( ^-f 1 + 2 Sn-lVl )(2a - 1)J . (5. 


16) 


The pressure distributions can be kept in their original form by replacing 
into them the expression ( i - 77 2 + 2it, 77 ) instead of ( i - 7? 2 ). The damped 
liquid force is 


F = m9 2 y e*^ 
y 




+ 2iC 


2n- 


-i^n- 


,) 


(5. 17) 


The moment with 



CO 



( a )< 2 n -‘, 4 


is 
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M = 
z 


4m a lfit 

— o— e 

7 r J o 


E 

n = i 


f( 2n- 1) tanh 

(<2n-i) 

4 

1 2 2 

a nh j 

' h a 

7T “ COSh 

a 

(2n-l)7rh 

^2n-l h 

7T 

a 

a 


a L. a 

- i + 2 «2„-A n -l) 


(5. 18) 


SECTION VI. NUMERICAL EVALUATION AND CONCLUSIONS 
A. Liquid Motion and Mechanical Analogon 

The natural frequencies of the liquid in the container are well de- 
scribed by the formula ( 2. 9) . They are indirectly proportional to the square 
root of the length of the container. Table 1 exhibits the frequencies that come 
into play during excitation in y-direction and z-direction if one chooses for the 
length of the container aj = 250 feet and a 2 = 150 feet respectively. Since no roll 
excitation is involved in the structural analysis only these frequencies are used. 
Since the liquid height is rather small compared to the length and width of the 
container the hyperbolic tangent function plays an important role. 

The circular frequencies, of which thirty are presented in Table 1, ex- 
hibit rather small values. In y-direction they range from cu* = 0. 22 to cu 30 = 3. 48 
rad/sec and in z-direction from = 0. 43 to w 3 o = 4. 83 rad/sec. Of course, the 
fundamental frequency uq = 0. 22 rad/sec for hj/aj or c uj = 0. 28 rad/sec for 
hs/aj and coj = 0. 43 rad/sec for hi/a 2 or cjj = 0. 53 rad/ sec for h 3 /a 2 play the 
most important role since the largest sloshing mass is connected with those. 

The values of the liquid height were chosen h^ = 10 feet, h 2 = 12 feet and 
h 3 = 16 feet. The sloshing liquid masses can be obtained from the mechanical 
analogon and are represented in Table 2 for various liquid heights and container 
excitation directions. The results show that the first sloshing mode is the most 
dominant and that eighty percent of the liquid mass sloshes in the first vibration 
mode. This is an extremely large value and requires special precautions. If 
four million gallons of water will oscillate, the system will experience very large 
forces and moments produced by the oscillating liquid. The sloshing mass of 
the second mode is only about eight percent of the total liquid mass and only one- 
tenth of that of the first mode, while the third sloshing mode exhibits only one to 
three percent of the total liquid mass participating in an oscillation. These mag- 
nitudes suggest, of course, that only the very large first sloshing mode has to 
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be considered for a stability and response analysis. The water in the container 
can, therefore, be described by a simple mechanical model consisting of a mass 
point attached with springs at the tank wall, such that the ratio of spring stiffness 
to sloshing mass k/m represents the square of the first circular natural frequen- 
cy oni of the liquid in the container. 

Before proceeding to the interaction of structure and liquid a few results 
of the fluid response to translatory excitation of the container are presented. 

The pressure distribution along the container bottom is given for various fluid 
heights and forcing frequencies in y- and z-direction. Figures 8 through ii 

exhibit the pressure (p, - Pgh)/pay e iSit . The maximum pressure is 

bottom o 

exhibited at left and right container walls and assumes for an excitation ampli- 
tude of y feet a value of about ( 55y + 62. 4) lbs/ft 1 2 , which yields for an exci- 
o o 

tation amplitude of y^ = i ft a value of 679 lbs/ft 2 . This happens if the forcing 

frequency is £2 = i. i oo * . For a larger fluid height of h 2 = 13 feet the maximum 
pressure at the container bottom for one foot excitation amplitude is about 
880 lbs/ft 2 * * . For excitation in z-direction where the container exhibits a length 
of a^ = 150 feet, the maximum pressure at the container bottom is 744 lbs/ft 2 
for a liquid height of hj = 10 feet and 960 lbs/ft? for a liquid height of h 2 = 13 feet 
(Fig. 11). Again, these pressures are based on one-foot excitation amplitude 
and a forcing frequency close to resonance, £2 = 1. loq. The response of the 
fluid force (lbs/ft) and liquid moment [ft - lbs/ft] are exhibited in Figures 12 
and 13. At resonance, of course, they exhibit singularities since no damping 
has been included in the results of the theory. For a very small damping a peak 
value of the force in y-direction of about 3. 63 x 10 8 lbs is obtained in the first 
resonance for a one-foot excitation amplitude, while the liquid moment exhibits 
at this location a value of about 6 x 10 8 ft-lbs about the center of gravity of the 
undisturbed liquid. 

B. Interaction of Structure and Liquid 

Here the free and forced oscillations have been treated. Further- 
more, the response of the system to various wind pulses has been performed. 

1. Free Oscillations. The knowledge of the coupled frequencies and 

mode shapes is of fundamental interest. In addition, the effect of the change of 

stiffness has been investigated and it was found that an increased stiffness would 

not pay off. Table 3 exhibits the values for the coupled natural frequencies 

coi; 2 of the sound suppressor system for the orientation in the z-direction. The 
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value a is a stiffness parameter exhibiting for unity in its magnitude the nominal 
stiffness and for a = 2 twice the stiffness of the support structure. 

It was found that the natural frequencies in y-direction for a liquid height 
of sixteen feet are extremely low ( 0.22 rad/sec) and that an increase in stiffness 
does not particularly pay off. In z-direction, where the sound suppressor system 
exhibits its largest area against a wind, the lowest natural frequency is presen- 
ted in Table 3 and is very close to the uncoupled fundamental natural frequency 
of the liquid. It is reached monotonically (w* = 0. 5319) with increasing stiffness 
of the structure. The second coupled natural frequency exhibits larger values 
(a little more than one cps for nominal stiffness of the support structure) and 
shows that a doubling of the stiffness of the support structure increases this 
frequency only by about forty percent. This frequency approaches with increas- 
ing stiffness the value infinity. 

2. Sinusoidally Forced Oscillations . Since a wind in the z-direction 
is more likely to build up in a regular fashion and since the area obstructing the 
wind perpendicular to the z-direction is the largest, this direction is chosen for 
the numerical evaluation of the response of the system to a sinusoidal excitation 
force. Although a pure sinusoidal wind of amplitude 13 m/sec is quite unlikely 
for a longer time period, the results will reveal some interesting facts about the 
influence of the damping of the structure and liquid. All wind data are based on 
the American Standards Association (ASA) paper A58. 1 - 1955, "American 
Standard Building Code Requirements for Minimum Design Loads in Buildings 
and other Structures," and the American Society of Civil Engineers (ASCE) 
paper No. 3269, "Wind Forces on Structures. " 

The response curves for which graphs have been omitted here show that 
the response of structure as well as liquid is very small as long as the forcing 
frequency O is not in the immediate neighborhood of one of the natural frequen- 
cies Wj or w 2 . An increase in structural damping decreases the response around 
the resonance peak of the structure considerably, but has hardly any noticeable 
influence upon the magnitude of the liquid response. The same is true for a 
change of liquid damping. An increase of the damping of the liquid will change 
(decrease) the peak length of the liquid response in the immediate vicinity of the 
natural frequency of the liquid, but does not effect the response of the structure. 

The following values exhibit the effect of the change of liquid and struc- 
tural damping upon the peak values for the deflection of the structure (zj max) 
and the amplitude of the liquid ( z 2 max) at the container wall. 
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£ = 0.001 

£ = 0. 01 

£ = 0.02 

£ = 0.05 


g=0. 001 g= 0. 02 

g=0.00i g = 0. 02 

g = 0.00i g= 0. 02 

g=0.00i g=0.02 

Z 1 

1 max 

50" 

2.5" • 

49. 5" 

2.5" 

49. 5" 

2. 5" 

49. 5" 

2. 5" 

z 2 

max 

9. 6" 

9. 6" 

5.3" 

5" 

3" 

3" 

1.25" 

1.25" 


Most likely the system will exhibit a liquid damping factor £ between £ = 0. 01 
and £ = 0. 03 and a structural damping of at least g = 0. 02. A maximum displace- 
ment of the container in the most unfavorable sinusoidal wind excitation case 
would be of a magnitude of 2. 5 inches, while the liquid height would exhibit a 
magnitude of about five inches. 


3. Wind Response . The response of the sound suppressor system 
has been numerically evaluated for various wind inputs in z-direction. An 
average value for the wind pressure of 50 psf has been employed. Figures 14 
through 37 exhibit the response of an exponentially decaying pulse for various 
stiffness parameters a and various degrees of decay. The equivalent wind force 
is F = 1. 1 x 10 6 lbs. It is assumed here that the system at the time t = 0 

Z eq * -fit 

was at rest and that at this time a wind of equivalent force F e hits the 

z eq. 

sound suppressor system. For /3 = 0 the step function is obtained (i. e. , a wind 


force which at t = 0 suddenly assumes a value F and remains at this value). 

z eq. 

The figures exhibit the displacement of the structure (z t /F ) in ins/lbs and 

Z 6C[* 

the liquid amplitude ( z 2 /2. 5F ) in ins/lbs , where the factor 2. 5 stems from 

Z 6C[* 

the conversion of the displacement z 2 to the liquid amplitude at the container wall. 
The main results are presented in Table 4. 


It can be seen that for. the case of excitation by a step function (/3 = 0) , 
the structure performs oscillations and exhibits a maximum displacement of 
about 2. 76 inches for half the nominal stiffness (a = %) . The liquid amplitude 
for these cases of stiffness are 7. 2 inches for half the nominal stiffness, 3. 6 
inches for nominal stiffness, 1. 77 inches for twice the stiffness and 0.7 inches 
for five times the structural stiffness of the support structure. 
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For |8 = 0. 1 a slight decay of wind is experienced. Here, the response 
of structure and liquid can be seen for the various stiffness cases in Figures 22 
through 29. The maximum structural amplitudes are 2 . 7 inches for half the 
nominal stiffness of the support structure, 1. 3 inches for nominal stiffness, 

0 . 66 inches for twice the stiffness and 0 . 28 inches for five times the nominal 
stiffness. The liquid exhibits a maximum amplitude of 5. 5 inches for o> = ^, 

2. 8 inches for a = 1, 1. 4 inches for a = 2 and 0. 55 inches for a = 5. 

Similar results are obtained for an exponentially decaying wind with 
/3 = 0. 5. This wind already reaches half its initial value in about 1. 6 seconds 
in contrast to the case 0 = 0. 1 where it takes 7 seconds for the wind to exhibit 
half its initial magnitude. The results of the response of the container and the 
liquid are exhibited in Figures 30 through 37. The maximum amplitudes reached 
for the structure are 2. 3 inches for a = 5 , 1. 2 inches for cn = 1, 0. 63 inches 
for a = 2 and 0. 23 inches for o; = 5. The maximum liquid amplitudes are 2 . 8 
inches for a = 1. 4 inches for a. = 1, 0. 7 inches for a. = 2 and 0. 28 inches 

for a = 5. The values in the parenthesis behind the maximum values in the 
table indicate the values which are assumed after the pulse has ceased to be 
effective. 

The response of the support structure and the liquid to a rectangular 
pulse of t = 2, 4 and 10 seconds duration is exhibited in Figures 38 through 61 
for various stiffness parameters (a = 1, 2 and 5). A wind is suddenly applied 

and ceases to act at a certain time t*. It may be noted that the pulse duration 
and the natural frequency have an important influence upon the magnitude of the 
vibration after the pulse has been completed. This can be seen from the follow- 
ing investigation. The solution for the motion of the system after the pulse has 
been completed is given by 


Zi = Ajlcos cu 1 (t - tj) - cos aqt] + A 2 [cos w 2 (t - tj) - cos co 2 t] 


z 2 = AjAjIcos u>j(t - tj) - cos uqt] + A 2 A 2 [cos co 2 (t - tj) - cos u> 2 t] 


where Aj and A 2 are some given amplitudes. Combining the terms in the paren- 
thesis to one trigonometric function yields 


cos oq(t - ti) - cos u^t = C sin (wjt - o;) 
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where 


co jtj co i 

C sin a = sin oofa = 2 sin — — cos 

lti 

C cos a = [1 - cos ojjtj] = 2 sin 2 


which gives the amplitude and phase relation 


C = 2 



CO it! 

and a = 


Therefore 


= 2Ai 

o>iti 

sin — 

sin [a; it - ] + 2A 2 

co 2 ti 
sm 2 

sin j cl> 2 t - w 2 t ! ] 



v 2 ' 


V 2 / 


z 2 = 

a) iti 

sm 2 

sin [wit - cuit : ]+ 2A. 2 A 2 

C0 2 tl 

sm 2 

sin { co 2 t - W 2 I 1 ) 



^ 2 ' 


V 2 ' 


From this it can be seen that the magnitude of the vibration after the pulse has 
been completed depends on the value j sin co jt x | and |sin u^tx j . If those values 

happen to be small values the response will be small. The main results are 
given in Table 4, while the more detailed response can be seen in Figures 38 
through 61. 

The response for a rectangular wind pulse of various durations 
(ti = 2, 4, 10 sec. ) for a = \ is exhibited in Figures 38 through 43. The maxi- 
mum amplitude of the container is 2. 76 inches in all three pulse cases. After 
the pulse has been completed the structure oscillates with about the same ampli- 
tude in the case of pulse duration ti = 2 seconds , while it drastically has changed 
its amplitude to a value of 0. 3 inches in the case of a pulse duration ti = 4 
seconds. For a pulse duration of t! = 10 seconds the amplitude remains 2. 76 
inches. The liquid amplitudes for these cases are in the same sequence given 
by 3. 6 inches, 6. 2 inches and 7. 2 inches. 
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After the pulse is completed the liquid amplitude remains about the same 
for the cases of pulse duration ti = 2 and 4 seconds, while for t* = 10 seconds it 
changes to the small value of 3. 5 inches. For nominal stiffness (a = 1) the 
results can be seen in Figures 44 through 49. For a pulse duration of tj = 2 
seconds the structure exhibits a maximum deflection of 1. 4 inches and assumes 
the value of 0. 6 inches after the pulse has been completed. For four seconds 
pulse duration the amplitude is 1. 4 inches; the structure then shows also a value 
of 1. 4 inches for the pulse duration of ten seconds, after which it performs 
oscillations of maximum amplitude of 1. 2 inches. The liquid exhibits in the 
first case of a pulse duration of 2 seconds the amplitude of 1. 8 inches , while for 
tj = 4 seconds, it is 3. 1 inches and finally for a pulse duration of tj = 10 seconds 
the amplitude of the liquid is 3. 6 inches during pulse duration and 1. 7 inches 
after the pulse has been completed. 

For twice the nominal stiffness (a = 2) the structure exhibits in all three 
pulse cases 0. 7 inches , while after the pulse has been completed the amplitude 
is 0. 01 inches for tj = 2 and 4 seconds and 0. 03 inches for t* = 10 seconds. The 
liquid exhibits values of 0. 9 inches throughout in the case t* = 2 seconds, 1. 5 
inches throughout in the case ti = 4 seconds, and 1. 8 inches during pulse dura- 
tion of length tj = 10 seconds and 0. 8 inches after the pulse has been completed. 
(See Figs. 50 through 55). 

For five times the nominal stiffness (a = 5) the structure exhibits in 
all three pulse cases the displacement 0.3 inches during the duration of the 
pulse, while after the pulse has been completed an amplitude of 0. 2 inches for 
ti = 2 seconds, 0. 3 inches for tj = 4 seconds and 0. 23 inches for t* = 10 seconds 
exists. The liquid exhibits an amplitude of 0. 36 inches throughout for t* = 2 
seconds, 0. 62 inches throughoutfor tj = 4 seconds, 0. 72 inches during pulse 
duration t* = 10 seconds, 0. 32 inches after the pulse has been completed 
(Figs. 56 through 61)'. 


Numerical Value? 
liquid mass m 


= 0. 865 x 10 5 


lbs sec 2 
inch 


liquid weight W = 3. 34 x 10‘ lbs 
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weight of support structure W = 5. 76 x io 6 lbs 

ss 

a lbs sec 2 

mass of support structure m = 1. 5 x 10 4 

ss inch 

weight of container W = 3. 58 x io 6 lbs 

c 

mass of container m = 0. 926 x 10 4 ~ ^ s sec 

c inch 

fluid heights: h* = 10 feet 
h 2 = 13 feet 
h'3 = 16 feet 

container dimensions: a* = 250 feet 

a 2 = 130 feet 

pairs of support beams: v = 16; p = 8 

Young's modulus of elasticity: E = 2. 1 x io 7 lbs/inch 2 

geometric moment of inertia of I = 7348 (inch) 4 

each pair of I-beams: f = lx io 6 (inch) 4 

z 

distance of the location of the 

sloshing mass from ground level: i = 90 feet = i. 08 x 10 s inches 

distance of the nonsloshing mass: L = 91 feet = 1. 09 x 10 3 inches 
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SECTION VI I. CONCLUSIONS AND RECOMMENDATIONS 


According to the previously presented results it can be seen that because 
eighty percent of the liquid (4 x 10® gallons) sloshes in the first vibration mode, 
the problem of interaction of structure and liquid can be a quite serious one. In 
y-direction the structure is a little weak and exhibits large deflections if excited 
in resonance. A regular pattern for a wind input is unlikely because the tower 
obstructs a clean build-up of winds; therefore, a slight stiffening is recommended. 
In z-direction a stiffening does not pay off, since the design of the support struc- 
ture seems adequate to withstand various wind pulses. 


George C. Marshall Space Flight Center, 

National Aeronautics and Space Administration, 

Huntsville, Alabama, September 30, 1'965. 
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natural circular frequency 


TABLE I. NATURAL CIRCULAR FREQUENCIES 



y- 

direction 

z - 

direction 


CO 

n 

hi/a t 

h 2 /a. t 


hi/as 

h 2 /a2 

h 3 /a 2 fluid height ratio 


“l 

0. 2248 

0.2258 

0.2832 

0. 4293 

0. 4683 

0. 5319 


co 2 

0. 4461 

0. 5051 

0. 5555 

0. 8354 

0. 9305 

1.0044 


co 3 

0. 6608 

0. 7422 

0. 8087 

1. 2024 

1.3104 

1. 3837 


“4 

0. 8663 

0. 9634 

1. 0380 

1. 5243 

1.6358 

1. 6851 


“5 

1. 0609 

i. 1666 

1. 2425 

1. 8028 

1. 8880 

1.9306 


“e 

1.2435 

1.3515 

1. 4236 

2. 0439 

2. 1104 

2. 1389 


CO 7 

1.4136 

1. 5191 

1. 5842 

2. 2548 

2. 3042 

2. 3223 


“8 

1. 5716 

1. 6710 

1.7276 

2. 4421 

2. 4774 

2.4886 


CO9 

1.7180 

1.8090 

1. 8569 

2. 6111 

2. 6358 

2. 6425 


“10 

1. 8538 

1.9353 

1. 9748 

2. 7660 

2.7829 

2.7869 


“11 

1. 9797 

2. 0515 

2. 0835 

2. 9098 

2.9213 

2. 9236 


“12 

2. 0970 

2. 1592 

2. 1848 

3. 0450 

3. 0526 

3. 0539 


“13 

2. 2066 

2.2598 

2. 2801 

3. 1730 

3. 1780 

3. 1788 


“14 

2. 3094 

2. 3545 

2. 3704 

3. 2951 

3. 2984 

3. 2988 


“15 

2. 4062 

2.4441 

2. 4565 

3. 4123 

3. 4144 

3. 4147 


“16 

2. 4979 

2. 5295 

2. 5390 

3. 5251 

3. 5265 

3. 5267 


“17 

2. 5851 

2. 6113 

2. 6186 

3. 6342 

3. 6351 

3. 6352 


“18 

2. 6683 

2.6899 

2. 6955 

3. 7400 

3.7406 

3. 7406 


“19 

2.7480 

2. 7657 

2. 7700 

3. 8427 

3. 8431 

3. 8431 


“20 

2.8247 

2. 8392 

2. 8424 

3. 9427 

3. 9429 

3. 9429 


“21 

2. 8987 

2. 9105 

2. 9129 

4. 0402 

4. 0403 

4. 0403 


“22 

2. 9703 

2. 9798 

2. 981 7^ 

4. 1353 

4. 1354 

4. 1354 


“23 

3.0379 

3. 0474 

3. 0488 

4. 2283 

4. 2283 

4. 2283 
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natural circular frequency 


TABLE I. (Concluded) 


y-direction z-direction 

hj/ aj h 2 / a t h 3 / a t hj/ a2 h 2 / a 2 h 3 / a 2 fluid height ratio 

w 2 4 3.1072 3. 1135 3. 1145 4.3192 4.3193 4.3193 

w 25 3.1730 3.1780 3. 1788 4.4083 4. 4084 4. 4084~’ 

w 2e 3.2372 3.2412 3. ¥418 4. 4956 4. 4957 _ 4. 4957 

u> 27 3.2999 3. 3031 T. 3036 4. 5813 4.~58i3 4.5813 

u> 28 3.3613 3.3639 3.3642 4.6654 4. 6654 _ 4. 6654 ~ 
w 29 3.4214 3.4235 374238" "477479 4. 7479 ~4. 7479 
w 30 3.4805 3.4821 3. 4823 4. 829l~4. 8291 4. 8291 ’ 
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TABLE H. SLOSH MASS RATIO 


n 

y- 

hi/a! 

•direction 

h 2 /a 1 

h 3 /ai 

z- 

hi/fy 

2n- 1/m 
-direction 
h 2 /a2 


i 

0. 806 

0. 803 

0.7998 

0.795 

0.785 

0.772 

2 

0. 086 

0.083 

0. 081 

0.077 

0. 070 

0.064 

3 

0. 029 

0. 027 

0. 025 

0.022 

0.019 

0.016 


TABLE IE. COUPLED NATURAL FREQUENCIES 
(a = STIFFNESS PARAMETER) 


a 

A z 

OJo 

Z z 

1 

2 

0. 528 

4. 76 

l 

0. 53 

6. 69 

2 

0. 53 

9. 43 

5 

0. 53 

14. 89 
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TABLE IV. WIND RESPONSE 


Decaying Pulse: 



z i 

z 2 

Z 1 

z 2 

z i 

z 2 

z i 

z 2 

0 

2.76" 

7.2" 

1.4" 

3. 6" 

0.7" 

1.77" 

0. 3" 

0.7" 

0. 1 

2.7" 

5. 5" 

1. 3"(0. 7") 

2. 8"( 1. 2") 

0. 66"(0. 38") 

1. 4"(0. 9") 

0. 28"(0. 14") 

; 0. 55"(0. 36") 

0. 5 

2. 3"( 1. 3") 









Rectangular Pulse: 


5 

1 

2 

X 

2 

1 

1 

2 

2 

5 

5 

z i 

z 2 

z i 

z 2 

Z 1 

z 2 

Z 1 

z 2 

2 


3. 6" 

1. 4"( 0. 6") 

1. 8" 

0. 7"(0. 01") 

0.9" 

0. 3"(0. 2") 

0. 36" 

4 

1 


6.2" 

1. 4"( 1") 

3. 1" 

. 0. 7"(0. 01") 

i 1. 5" 

0.3" 

10. 62" 

10 

2.76" 

7. 2"( 3. 5") 

1. 4"( 1. 2") 

3. 6"( 1. 7") 

0. 7"( 0. 03") 

1. 8"(0.8") 

! 0. 3"(0.23") 

0. 72"(0. 32") 























































•RDINATE SYSTEM 
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Exponentially decaying pulse 



t 



Exponential asymptotic step 










z x * structure in inches 
z 2 x 2.5= liquid in inches 
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FIGURE 15. RESPONSE OF LIQUID TO STEP FUNCTION (a = 
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FIGURE 16. RESPONSE OF STRUCTURE TO STEP FUNCTION (a = 1) 
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FIGURE 22. RESPONSE OF STRUCTURE TO EXPONENTIALLY DECAYING 
PULSE (/? = 0. 1) (a = i) 
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FIGURE 24. RESPONSE OF STRUCTURE TO EXPONENTIALLY DECAYING 
PULSE (a = 1) (0 = 0. 1) 
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FIGURE 26. RESPONSE OF STRUCTURE TO EXPONENTIALLY DECAYING 
PULSE (a =2) (/3 = 0. 1) 
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sec 


FIGURE 27. RESPONSE OF LIQUID TO EXPONENTIALLY DECAYING 
PULSE (a =2) (p = 0. i) 






FIGURE 28. RESPONSE OF STRUCTURE TO EXPONENTIALLY 
DECAYING PULSE (ce = 5) (0 = 0. 1) 
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FIGURE 31. RESPONSE OF LIQUID TO EXPONENTIALLY 
DECAYING PULSE ( a = |) (j8 = 0. 5) 
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FIGURE 32. RESPONSE OF STRUCTURE TO EXPONENTIALLY 
DECAYING PULSE (a = 1) (/3 = 0. 5) 
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FIGURE 35. RESPONSE OF LIQUID TO EXPONENTIALLY 
DECAYING PULSE (a = 2) (8 = 0. 5) 











sec 


FIGURE 37. RESPONSE OF LIQUID TO EXPONENTIALLY 
DECAYING PULSE (a = 5) <0 = 0. 5) 
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FIGURE 40. RESPONSE OF STRUCTURE TO RECTANGULAR 
PULSE OF DURATION U = 4 SECONDS (a = 1) 




sec 


FIGURE 41. RESPONSE OF LIQUID TO RECTANGULAR 

PULSE OF DURATION t x = 4 SECONDS (a = 1) 
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FIGURE 43. RESPONSE OF LIQUID TO RECTANGULAR 

PULSE OF DURATION tj = 10 SECONDS (a = A) 


100 









iiii i 1 1 1 1 1 1 1 1 1 1 u i 1 1 iii 1 1 1 1 j 1 1 i 1 1 J i ii | 

1 1 ! I T T J ! I ! ! ! ' I I I f I 




1 


t 




I I I I II- I -III HI III li- 

lt ■ ell i m h 1 1 n 1 1 1 n mi nr i ■ - I 

I { 1 1 I it I 1 I i i 1 i ‘ ill i i i 1 i i 

I I i I . J . l I ! L I I I I ! ( I | II j i : i I J i J 


ip * U 


III! 


I H i 




i/i i i i 






i j 1 1 1 1 1 1 1 1 1 1 klj/l j 1 1 j i in 1 1 1 1 n 


FIGURE 45. RESPONSE OF LIQUID TO RECTANGULAR 

PULSE OF DURATION (a = i) U = 2 SECONDS 
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FIGURE 46. RESPONSE OF STRUCTURE TO RECTANGULAR 
PULSE OF DURATION (a = 1) U = 4 SECONDS 
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FIGURE 47. RESPONSE OF LIQUID TO RECTANGULAR 

PULSE OF DURATION (a =1) tj = 4 SECONDS 


104 










Zx 



FIGURE 48. RESPONSE OF STRUCTURE TO RECTANGULAR 
PULSE OF DURATION (a = 1) t A = 10 SECONDS 
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FIGURE 52. RESPONSE OF STRUCTURE TO RECTANGULAR 
PULSE OF DURATION (a = 2) U = 4 SECONDS 
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FIGURE 53. RESPONSE OF LIQUID TO RECTANGULAR 

PULSE OF DURATION (a = 2) t t = 4 SECONDS 
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FIGURE 54. RESPONSE OF STRUCTURE TO RECTANGULAR 
PULSE OF DURATION (a =2) t t = 10 SECONDS 


111 




FIGURE 55. RESPONSE OF LIQUID TO RECTANGULAR 

PULSE OF DURATION {a. = 2) tj = 10 SECONDS 
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FIGURE 59. RESPONSE OF LIQUID TO RECTANGULAR 

PULSE OF DURATION (a = 5) ti = 4 SECONDS 
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FIGURE 60. RESPONSE OF STRUCTURE TO RECTANGULAR 
PULSE OF DURATION (a =5) = 10 SECONDS 
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"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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